Why do naked singularities form in gravitational collapse? 
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We investigate what are the key physical features that cause the development of a naked singularity, 
rather than a black hole, as the end-state of spherical gravitational collapse. We show that sufficiently 
strong shearing effects near the singularity delay the formation of the apparent horizon. This exposes 
the singularity to an external observer, in contrast to a black hole, which is hidden behind an event 
horizon due to the early formation of an apparent horizon. 
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I. INTRODUCTION 

In the past decade or so, several scenarios have been 
discovered where the gravitational collapse of a massive 
matter cloud results in the development of a naked sin- 
gularity iQ] . The final outcome of gravitational collapse 
in general relativity is an issue of great importance and 
interest from the perspective of black hole physics as well 
as its astrophysical implications. When there is a contin- 
ual collapse without any final equilibrium, either a black 
hole forms when the super-dense regions of matter are 
hidden from the outside observer within an event hori- 
zon of gravity, or a naked singularity results as the end 
product, depending on the nature of the initial data from 
which the collapse develops. 

The theoretical and observational properties of a naked 
singularity would be quite different from those of a black 
hole (see Q for further discussion of this). Thus it is 
of crucial importance to understand what are the key 
physical characteristics and dynamical features in col- 
lapse that give rise to a naked singularity, rather than a 
black hole. While many models of naked singularity for- 
mation within dynamically developing collapse scenarios 
have been found and analyzed not much attention 
has been given to understanding this important aspect. 
We begin here an investigation of this question. 

The main purpose of this paper is to identify the phys- 
ical process which exposes the singularity. We find that 
it is shearing effects which, if sufficiently strong near the 
central worldline of the collapsing cloud, would delay the 
formation of the apparent horizon so that the singular- 
ity becomes visible and communication from the very 
strong gravity regions to outside observers becomes pos- 
sible. When the shear is weak (and in the extreme case 
of no shear), the collapse necessarily ends in a black hole, 
because an early formation of the apparent horizon leads 
to the singularity being hidden behind an event horizon. 

For spherical gravitational collapse of a massive matter 
cloud, the interior metric in comoving coordinates is 



The matter shear is 



(2) 



where hab = gab+UaUb is the induced metric on 3-surfaces 
orthogonal to the fluid 4- velocity u", and is a unit 
radial vector. 

The initial data for collapse are the values on i = ij 
of the three metric functions, the density, the pressures, 
and the mass function that arises from integrating the 
Einstein equations (for details see e.g. jj]), 



(3) 



F{t;,r) = / p{t;,rydr, 



where 47ri^(ii, rb) = M, the total mass of the collapsing 
cloud, and where r > rb is a Schwarzschild spacetime. 
We use the rescaling freedom in r to set 



R{tu 



(4) 



so that the physical area radius R increases monotoni- 
cally in r, and with R[ — 1 there are no shell-crossings 
on the initial surface. (We will be interested here only 
in the central shell-focusing singularity at i? = 0, r = 
which is a gravitationally strong singularity, as opposed 
to the shell-crossing ones which are weak, and through 
which the spacetime may sometimes be extended.) The 
evolution of the density and radial pressure are given by 



F' 



Pr 
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The central singularity at r = 0, where density and cur- 
vature are infinite, is naked if there are outgoing non- 
spacelike geodesies which reach outside observers in the 
future and terminate at the singularity in the past. Out- 
going radial null geodesies of Eq. (|h are given by 
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dt 
dr 



(6) 



Consider first the case of homogeneous-density coi- 
lapse, p = p{t). Writing / = e^^^R'^ — 1, the Einstein 

-F/R. Then Eq. M) can 



equat 
be written 



ions give / 
itten as (5| 
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where m = r" (a > 1). If there are outgoing radial null 
geodesies terminating in the past at the singularity with a 
definite tangent, then at the singularity we have dR/du > 
0. For homogeneous density, the entire mass of the cloud 
collapses to the singularity simultaneously at the event 
{t ^U,r ^ 0), so that F/R oo. By Eq. (0), dR/du 
— oo, so that no radial null geodesies can emerge from the 
central singularity. It can be similarly shown that all the 
later epochs t > tg are similarly covered. 

We have thus shown that for spherical gravitational 
collapse with homogeneous density (and arbitrary pres- 
sures), the final outcome is necessarily a black hole. We 
note that this conclusion does not require homogene- 
ity of the pressures Pr and p±, and is independent of 
their behavior. The result generalizes the well-known 
Oppenheimer-Snyder result for the special case of dust, 
where the homogeneous cloud collapses to form a black 
hole always. 

An immediate consequence is that if the final outcome 
of spherical gravitational collapse is not a black hole, then 
the density must be inhomogeneous. In any physically 
realistic scenario, the density will be typically higher at 
the center, so that generically collapse is inhomogeneous. 



II. INHOMOGENEOUS DUST 

Consider now a collapsing inhomogeneous dust cloud 
(p = 0), with density higher at the center. The metric 
is Tolman-Bondi-Lemaitre, given by Eq. (0) with 
and e^'i' = i?'V(l + /), and 



i?2 = f{r) + 



F{r) 
R 



(8) 



These models are fully characterized by the initial data, 
specified on an initial surface t = ti from which the col- 
lapse develops, which consist of two free functions: the 
initial density pi(r) = p(ti,r) (or equivalently, the mass 
function F{r)), and /(r), which describes the initial ve- 
locities of collapsing matter shells. At the onset of col- 
lapse the spacetime is singularity- free, so that by Eq. (||), 



F{r) = r^F{r) , < F{0) < oo . 
The initial density pi{r) is 

Pi(r) = r-^F\r) . 



(9) 



(10) 



The shell-focusing singularity appears along the curve 
t = ts{r) defined by 



RMr),r) =0. 



(11) 



As the density grows without bound, trapped surfaces 
develop within the collapsing cloud. These can be traced 
explicitly via the outgoing null geodesies, and the equa- 
tion of the apparent horizon, t — iah(''), which marks the 
boundary of the trapped region, is given by 



i?(tah(r),r) =F(r). 



(12) 



If the apparent horizon starts developing earlier than the 
epoch of singularity formation, then the event horizon 
can fully cover the strong gravity regions including the 
final singularity, which will thus be hidden within a black 
hole. On the other hand, if trapped surfaces form suffi- 
ciently later during the evolution of collapse, then it is 
possible for the singularity to communicate with outside 
observers. 

For the sake of clarity, we consider marginally bound 
collapse, f — 0, although the conclusions can be gener- 
alized to hold for the general case. Then Eq. (^) can be 
integrated to give 



Ry^{t,r)^r^/^-l{t~-ti)F^/^ir) 
and Eqs. (|l|) and (ph lead to 



ts{r)=ti + 



F{r 
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iah(r)=ts(r)--F(r). 

The central singularity at r = appears at the time 

2 



ts(0) 



(13) 

(14) 
(15) 

(16) 



where where pc = Pi(0). Unlike the homogeneous dust 
case (Oppenheimer-Snyder), the collapse is not simul- 
taneous in comoving coordinates, and the singularity is 
described by a curve, the first point being (t = to,r = 0). 
For inhomogeneous dust, Eqs. (^) and ([l^) give 



^2 _ 1 ^ nb 



(3F - rF 



/n2 



GR-^R'^F 

A generic (inhomogeneous) mass profile has the form 



(17) 



F{r) = For^ + Fir^ + Far^ 



(18) 



near r — 0, where Fq — Pc/S- Homogeneous dust 
(Oppenheimer-Snyder) collapse has Fn = for n > 0, 
and Eq. (|l^) implies tr = 0. The converse is also true 
in this case: if we impose vanishing shear cr = 0, we 
get Fn — 0. Whenever there is a negative density gra- 
dient, e.g., when there is higher density at the center, 
then Fn ^ for some 7i > 0, and it follows from Eq. (oT 
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that the shear is then necessarily nonzero. Note that if we 
want the density profile to be analytic, we can set all odd 
terms F2„_i to zero; however, we note that this is not as 
such required by our own analysis, which is independent 
of any assumptions on Fn ■ 

The important question is: what is the effect of such 
a shear on the evolution and development of the trapped 
surfaces? In other words, we want to determine the be- 
havior of the apparent horizon in the vicinity of the cen- 
tral singularity at i? = 0, r = 0. To this end, let the first 
non-vanishing derivative of the density at r = be the 
Ti-th one (n > 0), i.e.. 

Fir) = For^ + Fnr''+^ + ■■■ , F„ < , (19) 

near the center. By Eqs. and @, 
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l-3F^^\t-U) + -Foit-Uf 

2n+l\ 



t^hir) = to - jFor'' 



F„. 



'iFc?'^ 



r" -f 0(r"+^). 



(20) 
(21) 



The time-dependent factor in square brackets on the right 
of Eq. (^0|) decreases monotonically from 1 dX t ~ io 
at t = tg. Thus the qualitative role of the shear in 
singularity formation can be seen by looking at the initial 
shear. The initial shear tJi = f7(ii,r) on the surface t = 
tj grows as r", n > 1, near r = 0. A dimensionless 
and covariant measure of the shear is the relative shear, 
|cr/9|, where 



e = 2— H — , 

R R' ' 

is the volume expansion. It follows that 



[1 + 0(r)] 



(22) 



(23) 



It is now possible to see how such an initial shear dis- 
tribution determines the growth and evolution of the 
trapped surfaces, as prescribed by the apparent hori- 
zon curve tahij"), given by Eq. (^l|). If we assume 
the initial density profile is smooth at the center, then 
Pi{r) = Pc + P2r^ + • ■ with p2 < 0, which corresponds 
to F{r) — Fqv^ + F2r^ + ■ ■ ■, with F2 < 0. Now suppose 
that p2 (and hence F2) is nonzero. Then Eq. (|l|) im- 
plies that the apparent horizon curve initiates at r = 
at the epoch Iq, and increases near r = with increas- 
ing r, moving to the future. Note that as soon as F2 
is nonzero, even with very small magnitude, the behav- 
ior of the apparent horizon changes qualitatively. Rather 
than going back into the past from the center, as would 
happen in the homogeneous case with F2 = 0, it is future 
pointed. This is what leads to a locally naked singularity. 
The singularity may be globally naked, i.e. visible to far- 
away observers, depending on the nature of the density 
function at large r. 



A naked singularity occurs when a comoving observer 
(at fixed r) does not encounter any trapped surfaces until 
the time of singularity formation, whereas for a black 
hole, trapped surfaces form before the singularity. Thus 
for a black hole, we require 



Uh{r) < to for r > 0, 



0. 



(24) 



In the general case (not necessarily smooth initial den- 
sity), this condition is violated for n — 1,2, as follows 
from Eq. (^). The apparent horizon curve initiates at 
the singularity r = at the epoch tQ, and increases with 
increasing r, moving to the future, i.e. iah > to for r > 
near the center. The behavior of the outgoing families of 
null geodesies has been analyzed in detail in these cases, 
and it is known that the geodesies terminate at the sin- 
gularity in the past [^|, which results in a naked singu- 
larity. In such cases the extreme strong gravity regions 
can communicate with outside observers. For the case 
n = 3, Eq. (24) shows that we can have a black hole if 



,5/ 



or a naked sing ularity, if F3 < -2F^ 
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^^3 > -2F,, , 
This is illustrated in Fig. 1. For n > 4, Eq. (|2^) is always 
satisfied, and a black hole forms. 
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FIG. 1. Apparent horizon curves near r = for the n — 3 
case, with _fo ~ 1- The labels on the curves give the values 
of ^3, the nonvanishing coefficient quantifying the shear. A 
black hole forms if > — 2. 

When the dust density is homogeneous, the apparent 
horizon starts developing earlier than the epoch of sin- 
gularity formation, which is then fully hidden within a 
black hole. There is no density gradient, and no shear. 
On the other hand, if a density gradient is present at the 
center, then the trapped surface development is delayed 
via shear, and, depending on the "strength" of the den- 
sity gradient/shear at the center, this may expose the 
singularity. It is the rate of decrease of shear as we ap- 
proach the center r = on the initial surface t — ti, 
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given by Eq. (^), that determines the end-state of col- 
lapse. When the shear falls rapidly to zero at the center, 
the result is necessarily a black hole; if shear falls more 
slowly, there is a naked singularity. It is thus seen that 
naked singularities are caused by the sufficiently strong 
shearing forces near the singularity, as generated by the 
inhomogeneities in density distribution of the collapsing 
configuration. When shear decays rapidly near the singu- 
larity, the situation is effectively like the shear-free (and 
homogeneous density) case, with a black hole end-state. 

It provides a useful insight to note that when a black 
hole forms, the apparent horizon typically springs into 
being as a finite-sized surface, at a finite r, then mov- 
ing to the center r = 0. This is what happens, for ex- 
ample, in the Oppenheimer-Snyder black hole formation 
in homogeneous dust collapse. In such cases, the event 
horizon, which does typically start at a point, could have 
formed earlier than the apparent horizon. On the other 
hand, in the case of a naked singularity, it follows from 
Eqs. and (|2^), that the apparent horizon starts at 
r — 0, and then is future directed in time, i.e. tah grows 
with increasing coordinate radius r along the apparent 
horizon curve R = F. These two behaviors of the ap- 
parent horizon curve are very different, and governed by 
shearing effects. A comoving observer will not encounter 
any trapped surfaces until the time of singularity forma- 
tion in the naked singularity case, whereas in the black 
hole case, the apparent horizon typically develops before 
the epoch of singularity formation. This is what we mean 
by delayed formation of the apparent horizon, caused by 
shearing effects. 

The relation between density gradients and shear may 
be understood via the nonlocal (or free) gravitational 
field. Density gradients act as a source for the electric 
Weyl tensor Q 

D'Eab^^BaP, (25) 

where is the covariant spatial derivative. (The mag- 
netic Weyl tensor vanishes for spherical symmetry.) In 
turn, the gravito-electric field is a source for shear (equiv- 
alently, the shear is a gravito-electric potential ^): 

2 2 

u'^VcCTafc + gQcTafc + CTac^r^'b " -^(^"^Kb = ~Eab ■ (26) 

Thus density gradients may be directly related to shear: 

Dap = -4c7Da(7 - 2eD^CT,t - 3D^ {u^Vc(Tab) 

- ioJ'Ti'abc " 3D^ {oaco^b) , (27) 

where we have used the shear constraint H^aab — fOa©- 
Equation (|27| ) makes explicit the link between the be- 
havior of density gradients and shear near the center, 
which was discussed above. The free gravitational field, 
which mediates this link, can also provide a covariant 
characterization of singularity formation. By Eqs. ( |2^ ) 
and ([2^), the relative gravito-electric field E/Q^ (where 
E"^ =~^E°-^Eab) near r = is given at t — ti by 



Thus naked singularities in spherical dust collapse are 
signalled by a less rapid fall-off of the relative gravito- 
electric field as we approach the singularity. Equa- 
tions (|2^) and (^8|) provide two equivalent ways of ex- 
pressing the result. This specifies how much shear is 
sufhcient to create a (locally) naked singularity. 

For the case of dust collapse, the role of shear in de- 
ciding the end-state of collapse is fairly transparent. To 
understand how shear affects the formation of the ap- 
parent horizon for general matter fields with pressures 
included is much more complicated, in particular since 
F = F{t,r), whereas F = for dust. In fact, even in 
some general classes of non-dust models (with nonzero 
pressure), it is possible to characterize collapse covari- 
antly. Above we showed that homogeneous density im- 
plies a black hole end-state. The next logical step would 
be to consider models for which the initial density is ho- 
mogeneous. For example, if the mass function is 

F{t, r) = f{r) - R\t, r) , /(r) = 2r^ , (29) 

then Eq. (||) shows that pi and (pr)i are constants. 
The density and pressure may however develop inho- 
mogeneities as the collapse proceeds, depending on the 
choice of the remaining functions, including in particular 
the initial velocities of the collapsing shells, and the col- 
lapse may then end up in cither a black hole or a naked 
singularity, depending on that (for a discussion on this for 
the case of dust collapse, we refer to [D). In fact, we can 
show that zero shear implies a black hole for these models. 
By Eqs. (||), (||) and (p9|), the shear- free condition leads 
to R' /R = l/r, and Eq. (||) then shows that p — p{t), 
i.e. the density evolution is necessarily homogeneous. As 
shown above, the collapse thus necessarily ends in a black 
hole. For the class of models given by Eq. (|29|), when- 
ever the collapse ends in a naked singularity, the shear 
must necessarily be nonvanishing. Although this class of 
models is somewhat special, the result indicates that the 
behavior of the shear remains a crucial factor even when 
pressures are nonvanishing. 

III. CONCLUSIONS 

Since black holes and naked singularities are of great 
interest in gravitation theory and astrophysics, it is im- 
portant to understand why these objects develop. The 
physics of this needs to be probed carefully in order to 
make further progress towards cosmic censorship, or to 
understand the physical implications of naked singulari- 
ties. 

It would appear that the only way a singularity can 
be laid bare is by distorting the apparent horizon surface 
and so delay trapped surface formation suitably. As we 
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have shown here, the shear provides a rather natural ex- 
planation for the occurrence of (locally) naked singulari- 
ties. Our main result is that sufficiently strong shearing 
effects in spherical collapsing dust delay the formation of 
the apparent horizon, thereby exposing the strong grav- 
ity regions to the outside world and leading to a (locally) 
naked singularity. When shear decays rapidly near the 
singularity, the situation is effectively like the shear-free 
case, with a black hole end-state. An important point 
is that naked singularities can develop in quite a natu- 
ral manner, very much within the standard framework of 
general relativity, governed by shearing effects. 

In the case of spherical dust collapse, shear and den- 
sity inhomogeneity are equivalent, i.e., the one implies 
the other. Although shear contributes positively to the 
focusing effect via the Raychaudhuri equation, 

&+\Q^ = ~\p-2a\ (30) 

its dynamical action can make the collapse incoherent 
and dispersive. (It is this feature which also plays the 
crucial role in avoidance of the big-bang singularity in 
singularity- free cosmological models 0|.) Depending on 
the rate of fall- off of shear near the singularity, its disper- 
sive effect can play the critical role of delaying formation 
of the apparent horizon, without directly hampering the 
process of collapse. The dispersive effect of shear always 
tends to delay formation of the apparent horizon, but 
is only able to expose the singularity when the shear is 
strong enough near the singularity. 

We have considered here spherical collapse. Very lit- 
tle is known about nonspherical collapse, either analyt- 
ically or numerically, towards determining the outcome 
in terms of black holes and naked singularities. How- 
ever, phenomena such as trapped surface formation and 
apparent horizon are independent of any spacetime sym- 
metries, and it is also clear that a naked singularity will 
not develop in general unless there is a suitable delay 
of the apparent horizon. This suggests that the shear 
will continue to be pivotal in determining the final fate 
of gravitational collapse, independently of any spacetime 
symmetries. In any case, our main purpose here has been 
to try to understand and find the physical mechanism 
which leads the collapse to the development of a naked 
singularity rather than a black hole in some of the well- 
known classes exhibiting such behavior. What we find is 
that the shear provides a covariant dynamical explana- 
tion of the phenomenon of naked singularity formation 
in spherical gravitational collapse. 
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